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Abstract 


The quest for the highest resolution microwave imaging and principle of time- 
domain imaging has been the primary motivation for recent developments in time- 
domain techniques. With the present technology fast time varying signals can now 
be measured and recorded both in magnitude and in-phase. It has also enhanced 
our ability to extract relevant details concerning the scattering object. In the past, 
the interface of object geometry or shape for scattered signals has received sub- 
stantial attention in radar technology. Various scattering theories were proposed 
to develop analytical solutions to this problem. Furthermore, the random inversion, 
frequency swept holography, and the synthetic radar imaging, have two things in 
common: (a) the physical optic far- field approximation and (b) the utilization of 
channels as an extra physical dimension, were also advanced. Despite the inherent 
vectorial nature of electromagnetic waves these scalar treatments have brought 
forth some promising results in practice with notable examples in subsurface and 
structure sounding. The development of time-domain techniques are studied 
through the theoretical aspects as well as experimental verification. The use of time- 
domain imaging for space robotic vision applications has been suggested. 
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I. Introduction 


The time-domain imaging is the synthesis of the scattered electromagnetic field 
distribution over an object plane. Over the past few years, an interest has developed 
in the use of time- domain methods or data in order to solve scattering problems. 
This approach has many attractive features, among them the large amount of in- 
formation contained in a single pulse response and simple physical interpretations; 
and improvement in time-domain technology increases the interest in such an ap- 
proach. The purpose of this paper is to present a short review of some techniques 
that have been used in the analysis of time-domain scattering. The problem of in- 
terest is the reconstruction of the shape of a convex scatter from knowledge of high 
frequency far field scattered from the object in response to a known incident field. 

The quest for the highest resolution microwave imaging and the principle of 
time-domain imaging have been the primary motivation for recent developments in 
time-domain techniques. In the last decade, modern techniques in sampling 
devices and the advent of fast pulse generators have brought new technology to 
the practicing radar engineers in the measurements of picoseconds at a greatly 
reduced cost. With the present new technology, fast time varying signals can now 
be measured and recorded both in magnitude and in-phase. 

The goals of NASA are to study the reflection coefficient at and below the sur- 
face of the tiles on the Space Shuttle with special emphasis on time-domain robotic 
vision applications and to continue research into problems related to the explora- 
tion of space and its applications of new technology and techniques. To obtain 
these goals, the Space Station shall be designed as a multi-purpose facility in which 
missions of long duration can be conducted and supported. These missions will 
include science and applications, observation, technology development and 
demonstration, commercial laboratories and production facilities, operational ac- 
tivities such as servicing/maintenance, repair of satellites, support of unmanned 
platforms, assembly of large space systems, and as a transportation node for trans- 
fer to other orbits and planetary missions. Automation and Robotic (A&R) is rapid- 
ly growing and it is an important technical area foreseen to increase productivity 
and enhance astronaut safety. The use of A&R for the Space Station can be viewed 
in two major areas: (1) teleoperated/robotic systems for servicing maintenance, 
repairs, and assembly, and (2) computerized systems to reduce the manpower re- 
quirements of planning, monitoring, diagnosis, control, and fault recovery of sys- 
tems/subsystems. 

In addition to increased productivity through autonomy, the A&R will result in 
increased operational capability and flexibility. Robotic operations for the Space 
Station will involve maintenance/repair of the entire structure including various sub- 
systems, orbiter/satellite servicing, astronaut assistance, equipment transfer, dock- 
ing and berthing, inspection, remote monitoring, rocket staging, telescience, and 
assembly of the station and large structures. To aid the astronaut in various tasks 
and replace him/her for some activities, robots must perform beyond the current 
state-of-the-art by responding to a high degree of environmental uncertainty and 
operational flexibility. In order to accommodate various performance goals in 
robots, design concepts have been proposed by Dr. Kumar Krishen of NASA and 
other researchers. 

Today’s technology enables us to extract pertinent details concerning the scat- 
tering object which conventional radar ranging lacks. The inference of object 
geometry or shape from scattered signals has received substantial attention in radar 
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technology largely due to its academic significant and the understandable commer- 
cial and military values. In the past decade, various scattering theories were 
proposed to develop analytical solutions to this problem. It is well known that ran- 
dom inversion, frequency swept holography, and the synthetic radar imaging have 
two things in common: (a) physical optic far-field approximation, and (b) utiliza- 
tion of the channels as an extra physical dimension. 

Despite the inherent vectorial nature of electromagnetic waves, these scalar 
treatments have brought forth some promising results in practice with notable ex- 
amples in subsurface and structure sounding. Recently, Dr. Kumar Krishen, of 
NASA/JSC, has proposed the use of time-domain imaging for space robotic vision 
applications. A multi-sensor approach to vision has been shown to have several 
advantages over the video approach. 

II. Theory 

The theoretical approach outlines the basic principles of time-domain imaging 
relative to the synthesis of the scattered electromagnetic field distribution over an 
object plane. It was stressed that it is the scattered field, not the total field that is 
of interest, it is known by virtue of the boundary conditions, the tangential electrical 
field must be zero at a perfect conducting surface, where as that of the magnetic 
field is discontinuous by the surface current. The literature review reveals that the 
theory will not change and it is applicable to any arbitrary source waveform. 

The scattering problems from rough surfaces is inherently different in nature 
from that of scattering by other bodies. Usually in the rough surface problem, an 
exact knowledge of the shape of the surface is neither available, nor is it of interest 
to the radar operator. Instead, only average properties of the surface shape enter 
into the problem. It is clearly understood that the last requirement rules out a bound- 
ary-value problem, since the exact boundary is not known. We are interested in 
the relationship between the average scattered field or radar cross section and the 
average surface properties. The radar cross section normalized by the area A, 
defined as, 

7(*<) = 7(*.M) (1) 

6i, 6s, and $s are shown for the bistatic case in figure 1 below. The rough surface 
is assumed to consist of height variations about a mean plane, which is taken as 
the xy-plane. The incident wave line in the sz-plane at a polar angle. 

Z 



Figure 1 - THE "BISTATIC SCATTERING GEOMETRY" 
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Figure 2- THE "SCATTERING GEOMETRY" 


Figure 1 shows the bistatic scattering for a planar rough surface whose mean 
height coincides with the xy-plane and figure 2 shows the scattering geometry. I is 
the plane of incidence and S is the scattering plane. 

The symbol < > refers to an average and the normalized average backscat- 
tering cross section are defined as 

7(*i) = 7(M,*) (2) 

Usually one has a specific purpose or reason for investigating scattering from 
a rough surface, these goals may be divided into three categories: (1) The problem 
of direct scattering; one wishes to know the average properties of the scattered sig- 
nal or cross section when the surface properties of rough surface are known, and 
the scattering information is expressed in terms of the surface properties. (2) The 
problem of inverse scattering; one wishes to obtain statistical information about the 
rough surfaces from a knowledge of the average properties of the scattered field. 
This problem is more difficult in that there appear to be many classes of rough sur- 
faces producing the same average scattering cross section as a function of the bis- 
tatic scattering angle and wavelength. 

The radar cross section (scattered field) can be expressed as 4* times the 
power delivered per unit solid angle in the direction of the receiver divided by the 
power per unit area incident at the target. The factor 4* enters from the definition 
of a solid angle. Assuming for the moment that the propagation path between the 
target and the receiving system is lossless, then this power ratio may be expressed 
as 


4*r a 


(£' • £••) 
(E* • £<•) 




(3) 


Where E s . and H s are the scattered electric and magnetic fields, respectively 
and E* and H‘are the incident fields. The scattered field is defined to be the dif- 
ference between the total field (with the target present) and the incident field. This 
is summarized as 


E s s E t - E*. 


(4) 
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Finding the cross section (a) now becomes a problem in electromagnetic field 
theory. In order that the cross section be independent of r, it is desirable to let r in 
Eq. 3 become arbitrarily large (scalar definition, not a tensor function). 


<r as 4» lim r J 

r-*oo 


(£*•£**) 

(£*£*•) 


. .. W-H”) 

4 ^ r (57T1F)* 


(5) 


To compute the scattered field, one sometimes computes the current induced 
on the target and then treats the target current distribution in terms of an equivalent 
aperture distribution. Antennae are often compared to an isotropic antenna, that 
is, anantenna which radiates uniformly in all directions. 


In the problem of clutter, returns from terrain are not wanted. However, in many 
cases the presence of clutter is unavoidable along with the desired signal. One can 
detect and analyze the desired signal significantly better if more is known about the 
properties of the clutter or noise produced by terrain scattering. It is noted that sur- 
face information is generally known, and the properties of the scattered signal are 
related to Eq. 3. 


The Fourier transform was used in this analysis and it was shown that space- 
time equations can be readily derived from their space-frequency counterparts. 
The equations derived will be applicable to field of arbitrary time variation. 


The source signal is related to the fourier transform, F(w), by: 


(<*) F{y>) = £j{t)J m *dl ( 6 ) m = 


( 6 ) 


The upper case lettering are used for the transformed functions unless other- 
wise noted. The wave number k = w/c, where c is the speed of light in the medium, 
a and p are the direction cosines of the vector r projected onto the xy-plane with 
a =sinecose, p=sinesine. 

It can be assumed that the field distribution over the antenna apertures and ob- 
ject plane is space-time separable; the field strength of each radiating element ob- 
served at some point r is e(r,t) = e(r) f(t - r/c), the field strength, and f(t - r/c) is the 
sine wave traveling at the speed of light. 

Assume the scatter in figure 2 to be placed in a far field, such that the assump- 
tion of uniform impinging planewave over the object can be obtained by the applica- 
tion of the far field Kirchhoff-Huygens principle. 

The basic principle of time-domain imaging over an object plane has been for- 
mulated and discussed by many researchers. It was not until 1983 when Wolfgang- 
M Boerner developed a finite difference-time domain method (FD-TD) with a near 
field to far field transformation using field equivalence, that a method to treat realis- 
tic scattering problems effectively was developed. In this method, the scattering 
problem is analyzed in two steps by treating the relatively complex near field region 
and the relatively simple far field region separately. The method involves first the 
determination of equivalence electric and/or magnetic current tangential to a virtual 
surface surrounding the scatter of interest by using the FD-TD method for a given 
external illumination. The computed near-field equivalent currents are then trans- 
formed to derive the far-field scattering problem and the radar cross section. Since 
the FD-TD method can deal with dielectric, permeable, and in homogeneous 
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materials in a natural manner, it is possible to incorporate most of the physics of 
wave interaction with any complex scatter of interest 

The Kirchhoff-Huygens principle deals with scattering, but the process is not 
as sensitive in analyzing scattering as the method put forth by Wolfgang-M Boerner, 
et al. TheKirchhoff- Huygens method deals with far field scattering, whereas 
Wolfgang- M Boerner treats both far field and near field scattering. 

The theory involving finite difference-time domain (FD-TD) is not without limita- 
tion. This method is based upon the physical optics approximation which is valid 
for scatter whose dimension is significantly larger than the incident wave length and 
is independent of polarization. 

: 1 



i i 


Figure 3 - SCHEMATIC DIAGRAM OF THE TARGET IMAGING SYSTEM 


a fCO *00 

/(to) = -jwF(w)T(w) I I e(r, y)e*p[— - a,x - 0,y)]dzdy. 
AXCT « «/-eo J ->oo 

/(to) = ^-jtvF(w)T(w)E t {ka t ,kp t )exp[-jkr.]. 

Alter 0 


( 7 ) 

( 8 ) 


Where Et (/fat Apt) is the transmitting-antenna anisotropic directional filter and r 
is the distance between the transmitter aperture and an object plane taken as 
reference. 

The second application of the Kirchhoff-Huygens principle and fixed phase 
evaluation can be shown with little calculus that the received waveform, g(t), is 

g(t) = const r F(w)T(w)K(w)[E„(ka 0 , *fl.)40*w)* £,(*»„ kp t )E r (ka r , k0r)E.(ka., kp.) e - ikr ] e~> mt dw. ( 9 ) 

where the direction cosine, the constant and the distance are dependent on the im- 
aging geometry. 

For the monostatic case commonly encountered (see figure 4 below), the 
transmitter and receiver coincide. In particular, if the phase center of the object is 
on boresight which suggests that the target is closely tracked, Eqn. 9 can be con- 
siderably simplified to yield: 

g(t) = const J~ F{w)T{w)R{w)[E a (ka a , kp.) + Ow) J £,(0, 0)£ r (0, 0)£.(*a„ kp.)e- 7 > ir ^ m, dw. ^ ^ 
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The delineation of the objective response is exercised as follows: First, the time 
response in the absence of any scattering objects is recorded. 


(«) 

to 



X = 


co = 


ox' + 0x/ 
i/or a +/ 3 J 
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Figure 4 - IMAGING GEOMETRY 
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(/) x = x coa — 

y s x'sintf -j-y'costf, 


(y) x' = xcostf-i- ysintf, 
j' = — x sin 0 -r y cos 6. 


(h) = cos(^ -*• tf), x = |OP|cos(^-rtf) = |OP|(cos^cos^ — sintfsintf). 


(i) and coa ^= j^pj 

THE ROTATION TRANSFORMATION 
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The antenna’s coupling filter is removed when the signal is subtracted from 
other records with the scatter present. Multiple path scattering between the object 
and surrounding obstacles can be gaited out in time when sufficiently large space 
is given. By definition, 

E,(ka„k0,) = r° r° e.(r, y)ezp[2jk(a Q z + (t,y))dxdy. 0 1 ) 

J-ooJ-oo 

The system kernel, K(w), can be obtained once the coupling term has been 
eliminated. 


K(w) = co*<(ii») a / , (w)r(ti;)ie(t«»)£,(0,0)E r (0,0)f: # (0,0)e-^‘ p V 


( 12 ) 


The use of mathematical techniques, provides an opportunity for the backscat- 
tering waveform to be simplified. To give the system kernel, K(w), and to obtain the 
electrical field strength in frequency domain and by the projection theorem, the 
response obtained is equivalent to a slice in the two dimension fourier plane (figure 
4), which is a representation 



Figure 4 - ONE-DIMENSIONAL TARGET FREQUENCY RESPONSE 


of the one-dimensional target frequency response in the two- dimensional fourier 
plane at a particular viewing direction. Equation 9 becomes 


E(ka, k0) s 


C(w)e*>* r » 

*(») 


(13) 


The application of the inverse fourier transform of equation 13 yields the object 
impulse response 

(14) 


r = «(< - 2x' tin -) 

J— CO C 


which can be considered as a series of scans across the object surface where each 
is an integral of the object field which is orthogonal to the viewing direction. The 
spatial resolution is thus achieved by the utilization of the time channel and the min- 
imal resolvable dimension ax is related to the sampling integral T, namely 

cT (15) 

“ 2«m(0) 


25-10 



a different impulse response is obtained by the rotation of the object plane about 
the z-axis, which corresponds to another slice in the Fourier plane. Further rota- 
tion over a 180 degrees span will fill up the entire Fourier plane. 

The two-dimension object field distribution is recovered, namely the recnstruc- 
tion strength, e(x, y) 

e(.,y) = - A J R ^E(ka,k0)MyK Vk{os * /,9) ^ (16 ^ 

a smoothing filter A(ka,kp) is introduced to compensate for the amplification of high 
frequency component due to the factor k in field strength equation [Eq. 16]. It is 
not necessary to have A, for noise free data. For practical numerical evaluation, 
the finite sampling rate will impose an upper limit, R, to the wave number k. To 
eliminate abrupt cutoff, which would cause ringing in the resultant image, it is neces- 
sary to introduce equation 16, to yield a smoothed image, equation 17 


A(ka,k0) = Min&£ 

(tt) 

«(*,») = E{ka,k0),in{ji'W ik W ) dk 

and is the basis of our reconstruction field strength. 


(17) 

(18) 


III. Experimental 

The government, private, and public sector have always had a need for time- 
domain imaging. This area of research did not surface until early 1962. The new 
and old researchers in electromagnetic field theory and other scientist provided 
many theoretical ideals on direct and indirect scattering. It was observed that these 
researchers did not have an insight into the circuit design, arrangement of ap- 
paratus and development of the necessary experiments to verify and determine the 
applications of the existing theory. 

Back scattering from metallic objects was experimentally observed and the sig- 
nal processing technique for extracting the object response was performed in 1985 
by Yeung and Evans. It was not until this time when Yeung and Evans put the ex- 
isting theory coupled with experimental verification together and made a projection 
of the long over due applications. After carefully reviewing the literature and the ex- 
perimental verification provided by these two researchers, it was observed that ap- 
plications of time-domain imaging is applicable to many aspects of the private and 
public sector, but the immediate and long range beneficiaries are NASA and the 
military. 

Recently, Dr. Kumar Krishen of NASA/JSC has proposed the use of time- 
domain imaging for space robotic vision application. It has been shown that space 
objects are covered with heat shields, (dielectric materials), and microwave 
penetrate through these shields. A multi-sensor approach to this robotic vision ap- 
plication has been shown to have several advantages over a video approach. 

The authors of this paper will use the necessary apparatus and experimental 
arrangement of Yeung and Evans to demonstrate the viability of impulse imaging. 
Verifying the experimental results should give us confidence in the validity of the 
theoretical method over a broad range of angles and the enormous experimental 
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potential of time-domain techniques in radar- scattering measurements. After 
scrutinizing the theory and experimental results of Yeung and Evans, the authors 
of this paper will look into time-domain robotic vision application and the effects of 
pulse width, polarization, look angles, and phase images of robotic vision coupled 
with, conducting experiments to verify the application of the existing theory and ex- 
periment. 

The target size employed in these experiments are comparable to the source 
pulse width, and it is well known that electrical shape of objects at low frequencies 
could be considerably different from the physical shape. 

Finite source rise-time (hence limited bandwidth) and insufficient source power 
are the major difficulties for achieving the highest resolution and sharp images. The 
system kernel involves the differentiation of the source waveform, and thisadds un- 
certainties in the determination of low-frequency components which is usually 
manifested as a slowly varying baseline in the deconvolved result. 

IV. Applications 

In 1986, a summary of goals and objectives were recommended by the NASA’s 
Advanced Technology Advisory Committee (ATAC) [1]. The robotic applications 
set up by the Committee are: (a) Teleoperation of mobile remote manipulator with 
collision avoidance, (b) Mobile multiple-arm robot with dexterous manipulators to 
inspect and exchange orbital replaceable units, (c) Systems designed to be ser- 
viced, maintained, and repaired by robots. 

The current remote manipulator system (RMS) will be used to manipulate ini- 
tial station assembly in the Space Shuttle. Further versions of the RMS will be 
employed to assemble Space Station and satellite. 

NASA/JSC is currently developing an early robot, the extra- vehicular 
astronaut (EVA) retriever, to perform as an aid to an EVA/man maneuverable unit 
(MMU) astronaut and retrieve the astronaut if circumstances dictate such an opera- 
tion. 


In the development of the space vision systems cost effectiveness, speed, small 
size, lightweight, high reliability and flexibility, and ease of operation must be con- 
sidered. Microwave techniques of time-domain imaging is one of the tools to 
achieve the above purposes. 

Robotic technology has been widely utilized by NASA to explore space; 
however, one of the most challenging problems to deal with is robotic vision. Vision 
is needed in fixed locations to keep track of location of objects, to pick and place 
machines/systems, and to monitor performance of non-automatic systems. In- 
deed, vision provides intelligence and flexibility in locating automated and non- 
automated machines. 

As the robotics era dawns in space, vision will provide the key sensory data 
needed for multi-faceted intelligent operations. In general, the 3D scene/object 
description along with location, orientation, and motion parameters will be needed. 
Sensor complements may include both active and passive microwave and optical 
types with multi-function capability. The fusion of the information from these sen- 
sors to provide accurate parameters for robots provides by far the greatest chal- 
lenge in vision. Furthermore, the compression, storage, and transmission of the 
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information associated with multi-sensor capability require novel algorithms and 
hardware for efficient operation. 

The vision requirements for space robotics are characterized by environmen- 
tal factors and tasks the robot has to perform. The natural space environment con- 
sists of intense light and dark periods. At a normal Space Station altitude, the 
sunlight intensity will fluctuate above 60 minutes of extreme brightness and 30 
minutes of darkness. Furthermore, due to the absence of atmosphere, light is not 
diffused/scattered. The ubiquity of white surfaces intensifies the problem of relying 
on photometric data for object identification/discrimination. A secondary source of 
concern effecting vision is the absence of gravity. For free flying and tethered ob- 
jects may be found due to the lack of disturbance caused by aerodynamic and 
gravitational forces. 

The unpredictable nature of maintenance and repair tasks creates a problem 
in the development of the capability/design of space robots. The vision capabilities 
must be adaptive/versatile to accommodate these uncertainties 

The quest for the highest resolution microwave imaging and the principle of 
time-domain imaging has been the primary motivation for recent developments in 
time-domain techniques. In the last decade, modern techniques in sampling 
devices and advent fast pulse generators has brought new technology to the prac- 
ticing radar engineers in the measurements of picoseconds at a greatly reduced 
cost. With the present technology, fast time varying signals can now be measured 
and recorded both in magnitude and in-phase. 

Conventional radar ranging has been enhanced by modern technology’s ability 
to extract details about the scattering object. 

V. Conclusions 

In this paper various research approaches to time-domain imaging have been 
reviewed with particular emphasis on the mathematical and physical approxima- 
tions made in them. 

Yeung and Evan dealt with the Kirchhoff-Huygens Principles, which deal with 
far field scattering whereas Boerner dealt with both far field and near field scatter- 
ing. We feel that the far and near field scattering technique is a better approxima- 
tion, particularly when applied to a scattering object that is not a perfect conductor. 

Based on the theoretical and experimental verification provided by these 
authors it is necessary that future research be done so as to go beyond the exist- 
ing approximations of Boerner, et al. For instance, a non-trivial generalization, 
should apply to conductivity and non-conducting targets. 
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